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Abstract 

We recall that the theory of electromagnetism consists of three building blocks: (a) the in- 
homogeneous Maxwell equations for the electric and magnetic excitations (D,H) (which reflects 
charge conservation), (b) the homogeneous Maxwell equations for the electric and magnetic field 
strengths (E, B) (which reflects flux conservation), and (c) the constitutive relation between (D, H) 
and (E,B). In the recent paper Lakhtakia proposed to change the standard boundary condi- 
tions in electrodynamics in order to exclude certain constitutive parameters. We show that this is 
inadmissible both from the macroscopic and the microscopic points of view. 
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I. INTRODUCTION 



Let us consider magnetoelectric matter with the constitutive relation 

D = eeoE + aB, (1) 
H = — B-aE. (2) 

Here Eq and /io are the electric and magnetic constants of the vacuum, e and /i are the 
permittivity and permeability of the matter, and a is the axion (Tellegen) parameter. 

In a recent discussion of the electrodynamics in magnetoelectric media of such a special 
type, Lakhtakia II] proposed to replace the standard boundary (jump) conditions on an 
interface between two media characterized by different values of the constitutive parameters 
e, fi, a by new boundary conditions which do not contain the axion parameter a. 

The motivation for such a replacement is as follows. When the medium is completely 
homogeneous, it is easy to verify that Maxwell's differential equations do not contain the pa- 
rameter a when (II])- ([21) is substituted into these equations. However, for the inhomogeneous 
situation, when there are two spatial domains with different values of a, for example, one 
should implement the boundary conditions on the surface that separates the two domains. 
The use of the standard boundary conditions allows for the influence of the nontrivial axion 
field a on the physical processes (in particular, on waves propagating through the interface). 
Lakhtakia l| modified the jump conditions such that this effect of a is removed. 

The authors of a comment j^] gave arguments against such an arbitrary modification 
of the boundary conditions. However, Lakhtakia 3|] disagreed, referring to the difference 
between the microscopic and macroscopic approaches to electrodynamics. Here we show 
that the claim of Lakhtakia is incorrect, both macroscopically and microscopically. 



II. MICROSCOPIC AND MACROSCOPIC ELECTRODYNAMICS 

We start from the microscopic Maxwell equations (see eqs. (l)-(4) of []J) that, after the 
spatial averaging, yield the equations of macroscopic electrodynamics: 

£oV-E = p, /io^ V X B -£oE = J, (3) 

V-B = 0, VxE + B = 0. (4) 
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Here we agree completely with Lakhtakia, namely that all four dynamical equations contain 
only two fields, E and B, and that all four equations hold in vacuum as well as in matter. 
There is however, an important difference: whereas in vacuum the electric charge and current 
densities p and J are defined by the free charges only, in matter the sources p and J are 
sums of both free and hound charges: 

p = pf + J = + J^ (5) 

The bound sources are related to the polarization and magnetization which provide an aver- 
aged description of the physical response of the microscopic constituents of matter exposed 
to the action of the electric and magnetic fields: 

pb = _v.p, j'^ = VxM + P. (6) 

Substituting ([5]) into ([3]), and making use of ([6]), we recast the inhomogeneous Maxwell 
equations into 

V • D = pf, V X H - D = (7) 
with the electric and magnetic excitation fields defined by 

D = £oE + P, H = po - M. (8) 

III. BOUNDARY CONDITIONS 

Let us consider the case when space is divided into two domains by a plane boundary 
surface S. When matter in each of the domains is homogeneous, it is convenient technically 
to study the electrodynamical processes (waves, in particular) in the two half-spaces sepa- 
rately. However, if we want the physical picture in the whole space, we will need boundary 
(jump) conditions for the electric and magnetic fields across the surface S. 

We derive here the boundary conditions in a slightly different way than it is done in 
textbooks. Namely, we will do it directly for the fundamental fields E and B and not for 
the excitations. Since each of the fields is a vector, we expect three boundary conditions for 
each field (one for the normal projection and two for the tangential components). Indeed, 
with the help of the usual technique by integrating the Maxwell equations ([3])-(jl]) in a thin 
region in a small vicinity of the surface iS, we derive the conditions for the fields at the 
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boundary: 

Ei^) - E« = .oV, E(^)-EW = 0, (9) 
B(2)-B«=M., B(2)-B« = 0. (10) 

The notation is obvious: the subscripts „ and r denote the normal and tangential projections, 
whereas the superscripts and ^'^^ label the half-space domains. 

The crucial feature of the boundary conditions across 5" is the presence of the surface 
charge and current densities Ps and in and (ITUI) . These surface densities have the form 

p, = pf+p^, J, = J: + J^ (11) 

thus including on an equal footing both the surface density of the free sources pi and 
and the surface density of the bound sources p^ and J^. The physical origins of the two 
types of the surface charges and currents on S are somewhat different. Whereas pi and 
describe the possible presence of the free sources right at the boundary (prepared under the 
conditions of an experiment, for example), the surface sources p^ and arise from the fact 
that matter (from the microscopic point of view) has different electromagnetic properties in 
the two half-spaces. As a result, the polarization and magnetization, while being continuous 
in each separate homogeneous domain, are discontinuous across the boundary S. 

We can easily derive the relation between the bound surface sources on 5* and the jump 
of the polarization and magnetization across S. Integrating the equations ([H]) in a infinitely 
thin neighbourhood including 5* (in the same way like for the Maxwell equations), we derive 

p[f^ - pi^) = -p's, M(2) - = J^ (12) 

Substituting ffTTj) and f|T2|) into fl9|)- f|T0|) . and taking into account ([8]), we recast the boundary 
conditions in a more familiar form: 

Di^) - D« = H(^) - hW = J^. (13) 

IV. AXION (TELLEGEN) MAGNETOELECTRIC MEDIUM AND THE POST 
CONSTRAINT 

Up to this point, we were quite general. Now, let us specialize to the case of the mag- 
netoelectric matter with the constitutive relation ([I])-([2])- The permittivity, permeability, 
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and the axion (Tellegen) field in the two half-space domains are £i,/ii,ai and 62,1^.2,0:2, 
respectively. Using ([I])- ([2]) in f|T3|) and fl9l)- f|T0l) . we find the boundary conditions 

Bo (e^Ei^) - eiE«) + («2 - «0B« = p^, (14) 
/^o ' {p2'bP - /^r^B«) - ia2 - aOE« = . (15) 

As we see, the axion field shows up explicitly. The analysis of the wave propagation in 
such a magnetoelectric medium reveals the influence of the axion field on the reflected and 
transmitted wave Q, 3 (see also 0,3). One can use this effect to measure the axion 
(Tellegen) field fl. 

Lakhtakia's proposal to change the standard boundary conditions in order to eliminate 



the contribution of axion (thus proving the so called Post constraint |4|) is physically unsub- 
stantiated. One cannot treat the boundary conditions as a kind of supplementary conditions 
that one can choose arbitrarily (like the the fixing of the gauge, for example). The boundary 
conditions (l9l)-( fT0|) and ( fT3l) are the Maxwell equations written in a different (integral and 
not differential) form for a specific region of space (i.e., for the infinitely thin neighbourhood 
of the boundary surface S). By changing a boundary condition, Lakhtakia actually changes 
the physical laws, namely, the Maxwell equations at the interface between the two media. 

An alternative (equivalent) explanation of the presence of the axionic terms in (fT^ - (fT5|) 
is as follows. Let us look at the right-hand sides of these equations. The free surface density 
sources pi and arise because, originally, in the Maxwell equations ([3]) there is a 5-function 
distribution of the free charge and current densities of the type = 5(x — x(S')) (and 
similarly for the free current density). The same applies also to the bound surface density 
sources. Indeed, the axion field, considered on the whole space, is a step function 

{«!, for X G 1st domain 
(16) 
a2, for X G 2nd domain 

Using this function in ([6]), we find a 5-function contribution from the magnetoelectric piece 
in the polarization and magnetization, namely, p^ = {ai — a2)Bn 5(x — x(S')) (and a similar 
expression for the bound current density). 

When these two delta-functions (one for the free and another for the bound sources) are 
substituted into the Maxwell equations ([3]), integration in the infinitesimally thin region 
around the boundary S yields the two contributions to the right-hand side of (!9|l- (|T0|) . The 



free source delta-function is responsible for the surface density terms p\ and J^, whereas the 
bound source delta-function gives rise exact 
argument was used in the previous comment 



to the axion terms in f[T4l) - f[T5l) . A similar 
|. In his response j^, Lakhtakia claimed that 



a microscopic approach and a "homogenization" of the fields might support his proposal. 
However, here we have analysed the problem starting from a microscopic viewpoint. It is 
unclear how any kind of "homogenization" can eliminate a delta-function at the boundary 
between the two domains filled with different matter. 

If we take any point at the boundary S and perform averaging and "homogenization" in an 
arbitrarily small neighourhood of this point, we will necessarily find two portions of space to 
the left and to the right of S", in which the electric and magnetic properties are homogeneous 
within the respective portions of the neighbourhood, but are not homogeneous and even not 
continuous across S. There just cannot be any "homogenization" across the boundary since 
S divides the two materials with essentially different physical properties. For example, 
we can have vacuum in the first half-space and a magnetoelectric medium in the second 
half-space. The vacuum is not polarized and magnetized. In contrast, the magnetoelectric 
medium becomes electrically polarized in a magnetic field and/or becomes magnetized in 
electric field, with the parameter a determining such polarization and magnetization. As a 
result, P and M are both discontinuous across S, and "homogenization" cannot change this 
fact. 



V. CONCLUSION 



Classical macroscopic electrodynamics (which can be consistently derived with the help 
of the spatial averaging from the microscopic electrodynamics) consists of three building 
blocks |8|: (a) the inhomogeneous Maxwell equations ([2]) (which reflects charge conserva- 
tion), (b) the homogeneous Maxwell equations (jlj) (which reflects flux conservation), and 
(c) the constitutive relation between the electromagnetic field excitations (D, H) and the 
electromagnetic field strength (E,B). The latter encodes the response of the medium to 
the action of the electric and magnetic fields in terms of the polarization and magnetization 
fields that are related to the bound charge and current source densities. Although some 
constitutive parameter (like the axion field a) may drop out of the differential equations ([3]), 
it still enters the constitutive law, reflecting the state of polarization of matter induced by 
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the magnetic field and/or the state of magnetization induced by the electric field. 

A modification of the boundary (jump) conditions across the surface 5* between the 
two different media, proposed in [II], is physically inadmissible because such a change of 
the boundary conditions amounts to a change of the Maxwell equations. Moreover, the 
averaging and "homogenization" arguments cannot eliminate the discontinuous behavior of 
the polarization and magnetization across the boundary S which is manifest in the delta- 
function like contributions both to the sources of the free charge and the bound charge. The 
Post constraint is unphysical and invalid. 
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